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We re-examine the classic problem of the rcnormalization of zero-point quantum fluctuations 
in a Friedmann-Robcrtson- Walker background. We discuss a rmmbcr of issues that arise when 
regularizing the theory with a momentum-space cutoff, and show explicitly how introducing non- 
covariant counter-terms allows to obtain covariant results for the renormalized vacuum energy- 
momentum tensor. We clarify some confusion in the literature concerning the equation of state of 
vacuum fluctuations. Further, we point out that the general structure of the effective action becomes 
richer if the theory contains a scalar field 4> with mass m smaller than the Hubble parameter H{t). 
Such an ultra-light particle cannot be integrated out completely to get the effective action. Apart 
from the volume term and the Einstein-Hilbert term, that are reabsorbed into renormalizations 
of the cosmological constant and Newton's constant, the effective action in general also has a term 
proportional to F{<p)R, for some function As a result, vacuum fluctuations of ultra-light scalar 

fields naturally lead to models where the dark energy density has the form pde(^) = pxit) + pz{t), 
where px is the component that accelerates the Hubble expansion at late times and pz{t) is an 
extra contribution proportional to H^{t). We perform a detailed comparison of such models with 
CMB, SNIa and BAO data. 



I. INTRODUCTION 



The computation of the expectation value of the 

energy-momentum tensor in curved space, and partic- 
ularly in a Friedmann-Robertson- Walker (FRW) back- 
ground, is a classic and much studied problem that has 
important implications for understanding the origin of 
dark energy. The rcnormalization of (0|T^^|0) in a FRW 
background can be performed either with covariant reg- 
ularization schemes, such as dimensional regularization 
or point-splitting [1], or using a cutoflF in momen- 
tum space [2 4]. The former option has the advan- 
tage of preserving general covariance explicitly, while a 
momentum-space cutoff can be more intuitive, for in- 
stance when developing naturalness arguments. The two 
schemes should give the same results for the renormal- 
ized quantities. However, the explicit breaking of general 
covariance by a momentum-space cutoff leads to some 
technical subtleties, that apparently have given rise to 
some confusion in both the old and the recent litera- 
ture. For instance it was claimed that regularizing a the- 
ory with a sharp momentum cutoff Ac in a FRW back- 
ground leads to pathologies [4] connected to the fact that 
the resulting vacuum energy density and pressure sat- 
isfy an equation of state (EOS) that is not consistent 
with general covariance [5]. It was also claimed that 
energy-momentum conservation is violated by the flow 
of momentum modes across the cutoff, induced by the 
cosmological red-shift [6-8]. There even appears to be 
some confusion as to what general covariance implies for 
(0]T^i^]0). For a constant term such as the A^ diver- 
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gence, general covariance implies that its contribution to 
(0|7)jjy|0) must be proportional to the metric g^,^, so its 
EOS parameter is w = —1. However, contrary to some 
statements in the literature (e.g. [8, 9]) this is not true for 
a time-dependent term such as the quadratic divergence 
H'^{t)Al, where H{t) is the Hubble parameter. We will 
show that general covariance implies that the contribu- 
tion of this term to (0]T^^]0) must be proportional to the 
Einstein tensor Gf^,^. The corresponding EOS parameter 
is time-dependent and equal to the total EOS parameter 
of the background, Wtot (i) = Ptot {t) / Ptot {t) , rather than 
being w = —1. 

The aim of this paper is two-fold. First, we will 
clarify the above technical issues. When using a non- 
covariant regularization scheme one must also include 
non-covariant counter-terms. We will see explicitly how 
introducing the appropriate non-covariant counter-terms 
allows us to obtain fully covariant results for the renor- 
malized quantities. As a consequence, no pathologies 
arise even when employing a sharp momentum cutoff. 
The renormalized energy density and pressure satisfy the 
standard energy-momentum conservation. The EOS pa- 
rameters associated to the terms arising from both the 
quartic and the quadratic divergences are dictated by 
general covariance and fully agree with the results ob- 
tained within covariant regularization schemes. We will 
also discuss a possible alternative to the standard rcnor- 
malization of vacuum energy, in which the subtraction of 
the flat-space contribution from the vacuum energy com- 
puted in FRW is performed already at the purely classi- 
cal level, invoking arguments from the ADM formalism 
of General Relativity, as well as from AdS/CFT corre- 
spondence, developing an argument that we presented in 
[10, 11]. 

After having shown that no new physics appear merely 

as a conseqiiencc of employing a non-covariant regulariza- 
tion scheme, we will turn to the second aim of this paper, 
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which is to show that, in contrast, the general structure 
of (0|T^iy|0) changes in an interesting manner if the the- 
ory contains a scalar field 4> whose mass m is smaller 
than the Hubble parameter H{t) at the time of interest. 
In particular, this applies to an ultra-light scalar field 
with m < Hq. In this case the vacuum expectation value 
(VEV) of the energy-momentum tensor T^^ acquires ex- 
tra terms, compared to the standard case where all mat- 
ter fields are very massive compared to H{t). In general 
the VEV of can be obtained by functional differen- 
tiation, 



(o|Vlo> 
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where S'off is the effective action for gravity, which is 
obtained by treating the metric g^v as a classical back- 
ground and integrating out the massive matter degrees of 
freedom (see e.g. refs. [12-15]). One can perform the cal- 
culation of the effective action using manifestly covariant 
regularization schemes such as dimensional rcgularization 
or point-splitting. The effective action is then explic- 
itly generally covariant and depends only on structures 
such as the volume term J d^Xs/^, the Einstein-Hilbert 
term J d'^Xy^—gR, and higher-derivative terms such as 
J d^x^—g R^i^R^'^ . The volume term is reabsorbed into 
a renormalization of the cosmological constant (giving 
rise to the cosmological constant naturalness problem) 
while the term in the effective action proportional to the 
Einstein-Hilbert term renormalizes Newton's constant. 
Higher-derivative terms such as R^i,R^^'^ can be seen as 
genuinely new effects due to vacuum fluctuations, but 
they can only be relevant when the curvature radius is of 
order of the Planck length. In particular, their impact is 
utterly negligible in the recent cosmological epoch. 

The basic idea we develop in this paper is that if a 
scalar field (j) with mass m < Hq is present it cannot be 
integrated out completely, since it is not massive with re- 
spect to the energy scale of interest in cosmology, H{t). 
Thus, after integrating over the massive degrees of free- 
dom, the eflfectivc action will be a functional of both the 
metric g^^ and the field <p (or, more precisely, of its low- 
frequency modes, that we will denote collectively as x). 
This allows for a richer structure of the effective action 
and therefore of the corresponding (0|T^j/|0). In partic- 
ular, the effective action will generically include a term 
oc J d'^Xy/—g f{x)R, with some function /(x)- Even at 
low curvatures such a term is not suppressed with re- 
spect to the Einstein-Hilbert term, and can have inter- 
esting phenomenological consequences. We will see that 
this indeed leads to an explicit realization of a model 
of interacting vacuum fluctuations recently proposed in 
refs. [10, 11], where the dark energy tracks the total en- 
ergy density of the Universe. 

The paper is organized as follows. In section II we dis- 
cuss the renormalization of vacuum energy in FRW using 
a momentum-space cutoff. In section III we discuss the 
effect of an ultra-light scalar field on (OlTpj^lO). We then 
construct a cosmological model based on these considera- 



tions and we analyze its phenomenological consequences. 

We use signature ( — h- 1--|-) and natural units where h = 
c = 1, so G = Mpj^, where Mpi is the Planck mass. When 
we specialize our results to FRW, we work in a spatially 
flat FRW metric with cosmic time t, scale factor a{t) and 
Hubble parameter H{t) = {l/a)da/dt = d/a. Today, the 
Hubble parameter and the critical density take the values 
Ho and po = 3HQ/{8nG), respectively. 



II. VACUUM ENERGY AND 
RENORMALIZATION IN FRW SPACE-TIME 

A. Renormalization and general covariance with a 
momentum-space cutoff 

To illustrate the kind of problems that arise when com- 
puting (0|Tj(i,y|0) with a cutoff Ag over momenta, let us 
begin with the divergence in flat space. Considering 
for instance a single massless real scalar field, a straight- 
forward calculation gives 



Pbare(Ac) = (0|Too|0) 



J4_ 
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(2) 



where the subscript "bare" emphasizes that this is still 

a bare, cutoff-dependent quantity, rather than the physi- 
cal, rcnormalized energy density. For the pressure, which 
is given hy p ~ (1/3) Xli(0|Tii|0): analogous compu- 
tation gives [4, 16] Pbare(Ac) = A^/(487r^), and therefore 



Pbare(Ac) = ^Pbare(Ac) . 



(3) 



This result has sometimes lead to the interpretation that 
the EOS of vacuum fluctuation is p = wp with w = 1/3 
(see e.g. footnote 19 of ref. [5]). In contrast, using a 
scheme such as dimensional regularization that preserves 
Lorentz invariance, in flat space one automatically gets 
(0|7)t!y|0) oc ry^jy and therefore p = —p. 

The resolution of this apparent discrepancy is that the 
energy density and pressure given in ecjs. (2) and (3) are 
bare quantities rather then rcnormalized ones and, in a 
regularization scheme that breaks Lorentz invariance, the 
EOS parameter Wbare = Pbare/Pbare is not the same as 
the physical EOS parameter for the rcnormalized quan- 
tities, Wren = Pron/Pron- Indeed, regularizing the theory 
with a cutoff over momenta breaks the Lorentz symme- 
try of Minkowski space, since the notion of maximum 
spatial momentum is not invariant under boosts. This 
means that in this scheme the renormalization procedure 
also involves counter-terms that are not Lorentz invari- 
ant. To give an explicit example of such a counter-term 
it is convenient to work with a generic metric g^^i, at 
first, before setting ^r^j, = r/^y at the end. This allows 
us to derive the vacuum energy-momentum tensor from 
(0|T^^|0) = -{2/y/^)SS/Sg'"', a procedure that facili- 
tates the extension to curved space-time. The standard 
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Lorentz-invariant counter-term in the action is propor- 
tional to the volume term 



Sa 



= ^(Ac) / 



(4) 



with ^(Ac) a suitably chosen function of the cutoff. An 
example of a Lorentz-breaking counter-term is 



(5) 



Using = —{'^/'2)^/^g^i,6g'^" we see that the 

contribution to (0|T^,y|0) from Sa is 



2 5Sa 



(6) 



while the contribution from Sb is 
2 SSb 00 



B{A,){g'^g,,-2Sl6'l). (7) 

Evaluating them on the flat metric = jy^j^ and us- 
ing the flat-space expression T^^^, = {p,p,p,p) (or, more 
generally, using the FRW metric g^i, = {—l,a^,a^,a^) 
and T^^ = (p, a^p, o^p, o^p)), we see that Sa produces 
a counter-term for the energy density given by pA = 
— ^(Ac) and a counter-term pA = —pA = ^(Ac) for 
the pressure. In contrast Sb produces the counter-terms 
Pb = +PB = -~B{Ac). 

For illustration we consider again a minimally coupled 
massless scalar field where pbarc(Ac) = A^/(167r^) and 
f bare (Ac) = A^/(487r^). After adding these countcr-terms 
to the bare A^ divergences we get 



Pren 



Pbaro(Ac) + Pcount(Ac) 

-|^-A(A,)-i3(A,), 

Pbare(Ac) +Pcount(Ac) 

^+^(A.)-5(Ac). 



(8) 
(9) 



We can now choose A{Ac} and B{Ac) such that A + B = 
A^/(167r^) — pfinite (where pflnite is a finite part) and A — 
B = -A4/(487r2)-pfi„ite, i.e. A{A^) = At/{48Tr^)~ p^^ite 
and S(Ac) = 2A^/(487r^). In this way we get p^en = 
Pfinito and pion = ^/Ofinitc, and therefore the rcnormalized 
quantities satisfy Prcn = — Prcn- Thus, Lorentz invariance 
can be restored for the renormalized quantities despite 
the fact that it was broken by the regularization scheme 
in the first place. 

Having understood this simple point in the flat-space 
case, we can apply a similar strategy in curved space, 
except that now the guiding symmetry principle is no 
longer Lorentz invariance, but rather general covariancc. 
In curved space one can again use regularization schemes, 
such as dimensional regularization or point-splitting, that 
preserve general covariancc explicitly, and in this case one 
finds a fully covariant result already at the level of the 



bare quantities. In particular, (0|T^j/|0) can be written 
as in eq. (1), where S'off- is made of coordinate- invariant 
quantities. For instance, for a minimally coupled scalar 
field in d = 4— e dimensions, the UV divergent part of the 
effective action has the form (see e.g. eqs. (6.44)-(6.49) 
of ref. [1]) 

5eff = J d^x^j[ci{e)+C2ie)R (10) 

+C3(e) {Rf.^paR^'''"' - R^.R^" - mR)] , 

where ci(e), C2(e), C3(e) contain a term diverging as 1/e 
and a finite part. It is interesting to compare this with 
the structure of divergences obtained using a cutoff Ac 
over momenta, and specializing for simplicity directly to 
a FRW background. In this case, considering for definite- 
ness a massless and minimally coupled real scalar field, 
we get ((0|T,^|0))bare = diag( 

Pbare j Pbare ) Pbare ) Pbare / ; 

where [4] 

Pbare(Ac) = ^ + ^^^^ + 0{H^ In Ac) , (11) 



Pbare (-^c 



167r2 167r2 

487r2 



■ ci- 



167r2 



0{H^ In Ac) .(12) 



The coefficient ci depends on the cosmological epoch: 
Ci = —1/3 during De Sitter (DS), Ci = -|-1 during ra- 
diation dominance (RD) and Ci = 2/3 during matter 
dominance (MD) [10]. 

The A^ divergence is of course the same as in flat space. 
Again, it cannot be renormalized simply using the covari- 
ant counter-term (4) , since the coefficient of the A^ term 
in the pressure is not equal to minus the coefficient of 
the A^ term in the energy density. Rather, one must 
also introduce the counter-term (5) and adjust these two 
counter-terms so that the corresponding renormalized en- 
ergy density and pressure satisfy Pren = — Pren- An equiv- 
alent way to phrase this point is to observe that the A^ 
terms in /9barc(Ac) and Pha.ie{Ac) do not satisfy the con- 
tinuity equation, i.e. 



Pbare + 3i? (pbare + Pbare) 7^ , 



(13) 



since for these quartically divergent terms pbarc = but 
Pbare + Pbare 7^ 0. This is just a conscqucnce of the 
fact that this regularization scheme breaks general covari- 
ancc. Therefore energy-momentum conservation, which 
is a consequence of general covariancc, does not hold for 
the bare quantities. There is nothing wrong or unphysical 
about this, as long as we choose the counter-terms such 
that the renormalized energy density and pressure satisfy 
energy-momentum conservation, which they do, since we 
saw that we can renormalize the quartic divergence so 

that, for this term, Pren = — Pren- 

Similar but somewhat more interesting considerations 
can be made for the quadratically divergent term. First 
of aU, let us see what the EOS of the term oc H^{t)A^ 
should be, in order to satisfy energy-momentum conser- 
vation. Writing p{t) = cH^{t), where c is a constant. 
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p{t) = w{t)cH'^{t), and plugging this into the conserva- 
tion equation p + 3H{p + p) = 0, we get 



3[1 + w{t)] 



1 



a dH^ 



(14) 



Using the Friedmann equation H'^{t) = 
{8TrG/3)J2iPi{t)j where the index i runs through 
all contributions to the energy density of the Uni- 
verse (matter, radiation, dark energy, etc.) and using 
Pi{t) (X a-3(i+"'*), we get 



w{t) = 



= Wtot{t) ■ 



(15) 



This is simply the total EOS parameter Wtot(i) = 
Ptot/ptot that tracks the EOS parameter of the dominant 
component, so it evolves from Wtot(i) — 1/3 during RD 
to Wtot(i) — during MD, and finally Wtot(i) — —1 if 
there is a dark energy dominated phase with wde — — 1- 
We can now see that the terms cx iJ^(t)A^ in cqs. (11) 
and (12) do not satisfy energy-momentum conservation, 
since the bare EOS parameter ci is not equal to u'tot(^)- 
Rather, the values ci = {-|-1, -1-2/3, —1/3} computed as- 
suming a pure RD, MD and DS epoch, respectively, are 
reproduced by 



ci(i) =Wtot(i)+2/3. 



(16) 



The situation is completely analogous to that discussed 
for the divergence. If the quadratically divergent 
terms in the bare energy density and pressure were satis- 
fying energy-momentum conservation we could rcnormal- 
ize them with a covariant counter-term, which in this case 
would simply be the Einstein-Hilbert term / d^Xs/^ R, 
so they would just be reabsorbed into a rcnormalization 
of Newton's constant. The fact that the quadratically di- 
vergent terms in the bare energy density and pressure do 
not satisfy energy-momentum conservation means that 
we must also introduce a non-covariant counter-term, 
chosen so that the contribution proportional to H'^{t) 
in the renormalized energy density and pressure satisfy 

Pron(i) = const. X H"^ {t) and Pren(<) = Wtot{t)p^enit). We 

explicitly construct the required non-covariant counter- 
term in appendix A. 

To conclude this section, let us observe that a dark 
energy term cx H'^{t)Mp^ can be obtained rather gen- 
erally, for instance from holographic considerations. In 
this context it is sometimes stated that, if the regulariza- 
tion of the divergent expressions is handled in a Lorentz- 
invariant manner, as with dimensional regularization, the 
corresponding VEV of the energy-momentum tensor will 
have the form (0|r^,^|0) oc gfj,„, i.e. w = —1 [9]. The dis- 
cussion presented above shows that this is not correct. 
First of all, Lorcntz invariance is irrelevant here as it is 
not a symmetry of the FRW space-time. Rather, the rel- 
evant symmetry is general covariance. As we have seen 
above, for the term oc H'^{t) in the energy density, gen- 
eral covariance fixes the EOS parameter to w = Wtot(i) 



rather than to w = —1. The two results coincide only in 
a De Sitter epoch. 

An equivalent way of understanding this point is to 

observe that a term oc H'^{t) in (0|T^i^|0) is obtained by 
taking the functional derivative of the Einstein-Hilbert 
term in the effective action (and not of the volume term) . 
Therefore its covariant form is proportional to the Ein- 
stein tensor G^^ rather than to the metric g^^, (see also 
the explicit covariant computation in [17]). Again, if 
H{t) is constant in time a g^j, and the two results 
agree, otherwise they are different. This is also the most 
direct way of deriving the EOS for the term oc H'^{t) 
in (0|r^i,|0): using the fact that this contribution to 
(0|Tpi,|0) is proportional to G^^, together with Einstein 
equations G^i, = (87rG)T*°' we sec that in the end this 
contribution to (0|T^^|0) is proportional to r*°*. So, in 
FRW, the term proportional to H^{t) in the vacuum en- 
ergy density and pressure satisfy Pvac/Pvac = Ptot/Ptot = 

Wtotit). 



B. Comparison with previous work 

The fact that the terms oc H'^{t)Al in Pbare(Ac) and 
Pbare(Ac) do uot rcspcct energy- momentum conservation 
appears to have generated some recurrent confusion in 
the literature. Rcnormalization of the energy-momentum 
tensor with a momentum-space cutoff in FRW was dis- 
cussed long ago in a classic paper by Fulling and Parker 
[4]. However, we believe that our discussion above sheds 
light on some aspects of ref. [4]. In that paper, the au- 
thors considered the (00) and the (Oi) components of the 
Einstein equations in a FRW space-time, both written in 
terms of the bare quantities^ 



:Ao 



An 



SttGo 
3 

-87rGoj3b 



Pbare , (17) 

(18) 



where Aq is the bare cosmological constant, not to be 

confused with the cutoff Ac, and Gq the bare Newton's 
constant. Regularizing the theory with a sharp cutoff in 
momentum space, and restricting for simplicity to the 
contribution of a minimally-coupled massless scalar field 
in De Sitter space, they find that both pbaro and pbare 
have a quartic as well as a quadratic divergence, given 
by eqs. (11) and (12), with ci = —1/3 for De Sitter. In- 
cluding also the finite parts pfinite and pfinite we therefore 



^ Compared to cq. (5.4) of rcf. [4], for notational simplicity we 
limit to the spatially flat case and we do not explicitly write the 
term which is related to the rcnormalization of terms quadratic 
in the Riemann tensor, denoted -ff^i/ in ref. [4]. 
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have 



Pbare(Ac) 
Pbare(Ac) 



167r2 167r2 



487r2 487r2 



+ 0(lnA,)+pfi„ite,(19) 
+ O(lnAe)+Pfi„ite.(20) 



Thus eq. (17) can be rewritten as 




Pfinite 



(21) 

which shows that, at least in this equation, the quartic 
divergence is reabsorbed into a renormahzation of the 
cosmological constant and the quadratic divergence into 
a renormahzation of Newton's constant. ^ 

The situation is however different in eq. (18) since, if 
we now insert the explicit expression of pbarej we see that 
in this equation neither the qiiartic nor the quadratic 
divergence can be reabsorbed into the above renormal- 
ization of Aq and Gq- The authors of ref. [4] propose the 
following solution. They observe that eq. (18) can be ob- 
tained by taking a time derivative of eq. (17), upon use of 
the conservation equation pbare + 3il (pbare + Pbare) = 0. 
However, as we saw above, Pbarc and Pbare given in 
eqs. (19) and (20), do not satisfy this conservation equa- 
tion, neither for the quartically divergent terms nor for 
the quadratically divergent terms. They then discard the 
results (19) and (20) arguing that a sharp cutoff in mo- 
mentum space is not a correct rcgularization, and that 
one should rather evaluate pbare using a smooth cutoff 
function, i.e. replacing 

/ dk^ dkf{k/aKc), (22) 
Jo Jo 

where f{k/aAc) is a smoothed version of the theta func- 
tion 6{l — k/aAc). (Here k denotes the comoving momen- 
tum, while in terms of the physical momentum, fcphys = 
k/a{t), the cutoff is at Ac.) Further, tliey argue that 

Pbare should be definedhy Pbarc + 3iJ(pbare+Pbare) = 0.^ 

In summary, according to ref. [4] the bare energy density 
should be computed as an integral over modes with a cut- 
off given by a smoothed theta function, rather than with 
a sharp cutoff, while the bare pressure, instead of being 
given by the corresponding integral over modes with a 



The logarithmically divergent terms, that we have not written ex- 
plicitly, are reabsorbed into the renormalization of higher deriva- 
tive operators such as R^^R^" [4]. 
3 Actually, they write ptare = I'Pbare = {'ia?a)~'^d{a^Pbs.ie)/d't, 
see cqs. (5.5) and (5.7) of ref. [4], which is the same as Pbare + 
3-ff(pbaro + Pbarc) = 0- The rolc of the smoothing function 
/(fe/aAc) is that the time derivative in the operator acts also 
on it, because of the factor a{t) in the argument of /, generating 
a term proportional to /', without which it would not be possible 
to impose eq. (13) to be satisfied. 



smoothed theta function, should be defined by the con- 
servation equation. 

Our discussion in the previous subsection however al- 
lows us to get a better understanding of this issue. 
First of all, there is nothing wrong with using a sharp 
momentum-space cutoff, and physical results should not 
depend on the choice of the regulator. The fact that 
it is impossible to reabsorb the divergences in eqs. (17) 
and (18) into a renormalization of Newton's constant (by 
means of a counter-term proportional to the Einstein- 
Hilbert term) is not a sign of a pathology. Instead, it is a 
direct consequence of the need for non-covariant counter- 
terms to restore the symmetry after regularizing with 
a scheme that breaks general covariance. The Fulling- 
Parker prescription, which consists of using eq. (22) and 

defining Pbare(Ac) from Pbare +3iJ(/C)bare + Pbarc) = 0, Can 

be seen as the definition of a regularization scheme that 
tries to enforce energy-momentum conservation by hand 
at the level of the bare quantities even in the presence 
of a momentum-space rcgularization. In this sense it is 
a legitimate, but certainly unnecessary (and somewhat 
akward) choice of regularization.^ 

A certain confusion seems to exist even in the more re- 
cent literature. When using a momentum-space cutoff it 
is not always appreciated that, even if the bare quantities 
do not satisfy energy-momentum conservation, 



Pbare + 3i/(pbare + Pbare) 7^ , 



(23) 



it is still possible to renormalize the theory such that the 
renormalized quantities satisfy it. 



3ff(pren +J'ren) = 0. 



(24) 



Equivalently, when the regularization scheme breaks gen- 
eral covariance one needs to introduce non-covariant 
counter-terms and as a consequence the bare EOS pa- 
rameter Wbare = Pbare/Pbare is different from the physical 
EOS parameter w = Pren/Pren- For instance, in ref. [7] 
the fact that the bare vacuum energy-momentum ten- 
sor does not satisfy the conservation equation leads the 
author to propose a modification of the Einstein equa- 
tions of the form G^v + E^,^ = STrGT^v, where Ej^i/ 
is a tensor that satisfies V|J,T,^"' = STrGV^T^"" ^ 0. 
However, once again, there is no special physical mean- 
ing in the fact that the bare energy-momentum tensor 



^ It should also be observed that general covariance is not fully 
recovered at the level of bare quantities even with this prescrip- 
tion. This clarifies another point of ref. [4] that the authors them- 
selves considered disturbing, namely the fact that some conver- 
gent terms in the expression for the bare energy density of order 
m^T~*, where m is the mass of the scalar field and T the param- 
eter of adiabatic regularization, do not form the (00) component 
of any divergenceless tensor, see the discussion at the bottom of 
page 186 of [4]. Once again, the point is that if the regulariza- 
tion scheme breaks general covariance, the bare quantities will 
not respect it cither. On the other hand, a generally covariant 
result for the renormalized quantities can be obtained upon use 
of non-invariant counter-terms. 
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is not conserved. What appears on the right-hand side 
of the Einstein equations reaUy is the conserved, physi- 
cal energy-momentum tensor constructed from the renor- 
mahzed vacuum energy density and pressure. 

Similarly, in ref. [8] the EOS of the term H'^{t)Kl 
was investigated and it was claimed that the EOS 
parameter for vacuum fluctuations is w = — 1 + 
A^/(97rM|[). However, this is a statement about the ra- 
tio Pbarc(Ac)/pbarc(Ac) and not about the ratio of the cor- 
responding renormalized quantities.^ Similar comments 
hold for the discussion of the EOS in [18]. 

The same criticism applies if one takes an effective 
field theory point of view rather than working in the 
context of full renormalization theory. Let us consider 
the theory which has been used to derive the result 
Pbaro(Ac) = A*/(167r^) as the low-energy limit of a more 
fundamental theory that takes over at a scale M. The 
total physical vacuum energy density will be given by 
the integral over the modes with k < M, that simply 
gives M^/(167r^), plus the contribution from the modes 
with k > M which can only be computed using the UV 
completion of the theory and will be of order M^, but 
comes with a numerical coefficient that cannot be com- 
puted with only the low-energy theory at hand. Thus, 
the exact coefficients of in the energy density and the 
pressure are unpredictable within an effective field theory 
and, therefore, the physical EOS parameter cannot be 
computed either. Thus, either one fixes the EOS param- 
eter by means of symmetry principles, i.e. by demanding 
Lorentz invariance for the A^ term or general covariance 
for the H^A^ term, or one simply has no handle on it. 
In any case, the EOS parameter derived from the ratio 
Pbare(Ac)/pbare(Ac) has no physical meaning. 

C. An alternative to cosmological constant 
renormalization: subtractions in classical GR 

The renormalization of the A^ divergence gives rise to 
the infamous cosmological constant naturalness problem: 
the divergence in the bare vacuum energy, pbare(Ac) = 
(0|Too|0), is canceled by a counter-term Pcount(A-c) that 
corresponds to a counter-term proportional to J d'^x 
in the action. The physical, renormalized vacuum energy 
is then p^-cn = /Obaro(Ac) + /3count(Ac). Thc naturalucss 
problem arises from the fact that Pcount(Ac) must cancel 
the A^ divergence in /Obare(Ac), leaving a finite part that 
is exceedingly small compared to A^, for all values of thc 
cutoff larger than, say, the TeV scale, where quantum 
field theory is well tested. 

At first sight the necessity of an extremely fine-tuned 
cancellation between pbare(Ac) and Pcount(Ac) can be dis- 



^ Futhcrmorc, thc derivation in ref. [8] assumed that the EOS pa- 
rameter consistent with general covariance is «; = — 1 while, as 
we saw above, for the term oc H^{t)A^ general covariance gives 

U) = Wtot{t). 



turbing. It should however be observed that neither the 
bare term nor the counter-term have any physical mean- 
ing and only their sum is physical, so this fine tuning 
is different from an implausible cancellation between ob- 
servable quantities. The same kind of cancellation ap- 
pears for instance in the Casimir effect, and this might 
provide a first hint of what could be the correct treat- 
ment of vacuum energies in cosmology. In the Casimir 
effect one computes the vacuum energy density of a field 
in a finite volume (e.g. the electromagnetic field between 
two parallel conducting plates at a distance L) and sub- 
tracts from it the vacuum energy density computed in 
an infinite volume. If we regularize with a cutoff Ac over 
momenta both terms diverge as A^, but their difference is 
finite and depends only on the macroscopic scale L. Thus 
the quantity of interest is the difference between the vac- 
uum energy in the given geometry and thc vacuum energy 
in a reference geometry, which is just flat space-time in 
an infinite volume. This might suggest that, to obtain 
thc physical effect of the vacuum energy density on the 
expansion of the Universe, one should analogously com- 
pute the vacuum energy density in a FRW space-time 
and subtract from it the value computed in a reference 
space-time, which could be naturally taken as Minkowski 
space. This procedure leads to a sort of "cosmological 
Casimir effect" . 

An immediate objection to this analogy could be that 
in special relativity the zero of the energy can be chosen 
arbitrarily and only energy differences with respect to the 
ground state are relevant. In contrast, in GR we cannot 
choose thc zero of thc energy arbitrarily. One typically 
expects that "every form of energy gravitates" , such that 
the contribution of Minkowski space cannot simply be 
dropped. While it is certainly true that in GR the choice 
of the zero of the energy is not arbitrary, thc point that 
we wish to make here is that the correct choice can be a 
non-trivial issue (see also [10, 11]). As a first example, 
consider the definition of energy in GR for asymptotically 
fiat metrics. To carefully define the energy associated 
with a given field configuration it is convenient to use the 
Hamiltonian formulation of GR, which goes back to the 
classic paper by Arnowitt, Dcscr and Misner (ADM) [19] 
(see e.g. the textbook [20] for a very clear recent review). 
In order to define the Hamiltonian of GR one must first 
work in a finite three-dimensional volume, where it takes 
the form 

^^GR = -ffbulk + -^boundary • (25) 

The bulk term i?buik is given by an integral over the 
three-dimensional finite spatial volume at fixed time, and 
-ffbouuciary by an integral over its two-dimensional bound- 
ary. If we try to define the energy of a classical field 
configuration as the value of this Hamiltonian evaluated 
on thc classical solution we encounter a problem: as 
a consequence of the invariance under diffeomorphisms, 
the volume term i?buik vanishes when evaluated on any 
classical solution of the equations of motion. Then the 
whole energy comes from the boundary term. However, 



7 



when the boundary term is evaluated on any asymptoti- 
cally flat metric (including Minkowski) it diverges as the 
boundary is taken to infinity. The solution proposed by 
ADM is to subtract the same boundary term evaluated 
on Minkowski space-time. The resulting energy (or mass) 
is finite and provides the standard definition of mass in 
GR, known as the ADM mass. For instance, when ap- 
plied to the Schwarzschild space-time, the ADM mass 
computed in this way turns out to be equal to the mass 
M that appears in the Schwarzschild metric. The ADM 
prescription can be summarized by saying that, in GR, 
the energy E associated to a classical asymptotically- 
flat space-time with metric g^v can be obtained from the 
Hamiltonian -^gr through 

E = HG^\g^,] - Hgr['?m.] , (26) 

where r/^j/ is the flat metric. Even if the context in which 
this formula is valid, namely asymptotically flat space- 
times, is different from the cosmological context that we 
are interested in here, eq. (26) still suffices to make the 
point that the intuitive understanding that GR requires 
any form of energy to act as a source for the gravitational 
field is not really correct. Equation (26) tells us that the 
energy associated to Minkowski space does not gravitate. 

The ADM prescription has also been generalized to 
space-times that are not asymptotically flat, subtract- 
ing the contribution of some reference space-time whose 
boundary has the same induced metric as the geometry 
under consideration [21-24]. The issue has been stud- 
ied in detail in particular for asymptotically AdS space- 
times, where the problem appears because the boundary 
term in the Hamiltonian diverges as the boundary used 
in its deflnition approaches the boundary of AdS. In this 
case it has been shown that the appropriate subtraction 
can be performed without even introducing a reference 
background, but by adding some local counter-terms to 
the boundary action. These are given by coordinate- 
invariant functional of the intrinsic boundary geometry, 
flxed by the requirement to obtain a flnite energy as the 
boundary is taken to infinity [25, 26]. This prescription is 
particularly appealing since for AdS the structure of the 
divergences of the boundary action is such that they can 
be removed by adding a finite polynomial in the bound- 
ary curvature and its derivatives. This is analogous to 
the fact that in QFT the UV divergences are removed by 
adding counter-terms that are polynomials of finite or- 
der in the fields. Indeed, in the context of the AdS/CFT 
correspondence, this way of removing divergences in the 
gravitational action on the AdS side corresponds to the 
renormalization of the UV divergences in the conformal 
quantimi field theory that lives on the boundary [25- 
29]. This subtraction procedure correctly reproduces the 
masses of various known space-times that are asymptot- 
ically AdS. For instance, in 2 + 1 dimensions it correctly 
gives the mass of the BTZ black hole (which is asymp- 
totically AdSs) and also reproduces the transformation 
law and conformal anomaly of the stress tensor in the 
dual CFT. In 3 -|- 1 dimensions it correctly reproduces 



the known value of the mass of the AdS4-Schwarzschild 
solution and similar results can be obtained for asymp- 
totically AdSs solutions. Furthermore, while this sub- 
traction procedure provides a zero mass for pure AdS4 
space-times, it gives a non-vanishing value for the mass 
of pure AdSs ; which exactly matches the Casimir energy 
of the dual J\f = A super Yang-Mills theory that lives on 
the global AdSs boundary with topology x R [25] . 

What wc learn from these examples is that the intu- 
itive argument that "GR requires that any form of energy 
contributes to the gravitational field" is in fact loose and 
generally incorrect. Rather, in all cases an appropriate 
subtraction, formulated either as the subtraction of the 
contribution of a reference space-time, or in terms of lo- 
cal counter-terms, is required already at the classical level 
in order to obtain an energy-momentum tensor which is 
well-defined and reproduces the known properties of the 
space-time under consideration. 

Prom this vantage point it is rather natural to assume 
that the definition of energy associated to zero-point 
quantum fluctuations in a curved background should not 
escape this general rule. As we have seen, regularizing 
the theory with a cutoff Ac over comoving momenta, in 
FRW the bare energy density of zero-point fluctuations 
takes the form 

[/Obare(Ac)]FRW = [pbaro(Ac)]Mink + C (-ff^(OAc) 

+C(ff*(t)lnAe) . (27) 

Here [/3baro(Ac)]Mink IS the bare vacuum energy density 
in Minkowski space and, for a field of mass m, its general 
UV structure is 

[Pbare(Ac)]Mink = ^(A^) + ©(m^A^) 

+0{m'^ In Ac) + finite part . (28) 

For example, for a real minimally coupled massless scalar 
field, [pbarc(Ac)]Mink = CiAj and 

[/Obare(Ac)]FRW = CiA^ + C2H^t)Al + O {H\t) Ink,) , 

(29) 

with ci = C2 = l/(167r^).® The quartic divergence is the 
same as that found in flat Minkowski space. This can be 
seen from the explicit computation and is of course more 
generally a consequence of the fact that the result must 
reduce to that obtained in flat space in the limit H ^ 0. 

Our discussion of the definition of energy in GR sug- 
gests the following alternative to the standard renor- 
malization of the cosmological constant. Rather than 



Actually, the value of C2 also depends on the choice of modes, 
and hence of vacuum state, for the field under consideration. The 
value C2 = l/(167r^) is obtained using the modes that correspond 
to the usual Bunch-Davies vacuum. The most general result is 

obtained by taking a combination of these modes with arbitrary 
Bogoliubov coefficients. This changes the numerical value of the 
constant C2, but not the general structure of divergences dis- 
played in eq. (27), see e.g. appendix B of [10]. 
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eliminating the divergence with standard renormal- 
ization, which allows us to get rid of it at the price 
of the naturalness problem, another option is that we 
should first perform the appropriate subtraction, corre- 
sponding to the correct choice of the zero of the energy 
density. To understand what is the correct subtraction in 
FRW, we consider the Friedmann equation when on the 
right hand side of the Einstein equations we use as source 
T^,u + (0|7;„.|0), where Tff = diag(-/3,p,p,p) is the ordi- 
nary classical contribution of matter, radiation, etc, and 
(0|T^,^|0) is the corresponding contribution of zero-point 
fluctuations. The Friedmann equation then reads 

i^'(^) = ^(p+(o|^oo|o)) . (30) 

We fix the subtraction procedure by requiring that 
Minkowski space (i.e. H{t) = 0) should be a solution in 
the limit p — >■ 0. This implies that all terms in [pvac]FRW 
that do not vanish when — >■ must be subtracted, 
i.e. we must subtract the vacuum energy computed in 
Minkowski space. This eliminates the divergence 
(as well as other flat-space divergences such as terms 
DC rn^A^ that appear for a massive field with mass m) 
in a way that docs not suffer from the naturalness prob- 
lem. The remaining divergences are renormalized in the 
standard way. In particular, the term oc H'^{t)K^ in the 
energy density is reabsorbed into the rcnormalization of 
Newton's constant (plus non-covariant counter-terms for 
renormalizing the pressure, as discussed in section II A). 
The rcnormalization of these remaining divergences does 
not suffer from any naturalness problem (and, in any 
case, they could not be affected by a constant shift in 
the zero-point of the energy).'' 

It is also interesting to observe that this procedure pro- 
vides a now perspective for interpreting the results of 
ref . [30] . There the authors show that in a theory with a 
large number M of quantum fields, because of the finite 
volume fluctuations of vacuum energy around its average 
value, Minkowski space would be unstable to black hole 
formation in regions of size 0{N^/H-p\), where l-p\ is the 
Planck length. This is of course a paradoxical result, and 
to escape this conclusion the authors of ref. [30] suggest 
that the quantum gravity length-scale for N fields should 
be of order M^^H-pi instead of simply l-p\. In other words 
the UV cutoff should be parametrically smaller than the 
Planck mass, and of order M~^/'^M-p\. If this was the 
case, semiclassical gravity would break down at this scale, 
and the calculation used to prove black hole formation 



Needless to say, none of the experiments that are usually men- 
tioned as probes of vacuum fluctuations, such as the Casimir 
effect or the Lamb shift, arc afTcctcd by such a choice of the zero 
point of the energy. These experiments, of course, do not probe 
directly the choice of the zero in the definition of the energy; 
rather they always probe energy differences, while the zero of 
the energy is defined, arbitrarily, by normal ordering the Hamil- 
tonian. 



would no longer be valid. With our proposal, in con- 
trast, the problem is eliminated simply because vacuum 
fluctuations in Minkowski space do not gravitate. 

By means of this approach it is also possible to un- 
derstand why phase transitions in the early Universe do 
not produce an exceedingly large vacuum energy den- 
sity [31]. Consider for instance a scalar field (f under- 
going spontaneous symmetry breaking (SSB), with the 
typical double-well potential 

v{^) = -\^J?^^ + \^\ (31) 

with fx^ > 0. Computing the energy density in the FRW 
background, one finds the vacuum energy associated to 
SSB to be 

PsSB = -J^-^(H + 2i72). (32) 

The first term is the fiat-space result that would typ- 
ically give an unacceptably high value for the vacuum 
energy density. For instance, the electroweak transition 
would give pggg = 0(100 GeV). With our proposal, how- 
ever, this Minkowski space term does not contribute to 
the gravitating energy density, and one remains with the 
second term in eq. (32), that is 0{iJ,^H^/X) and therefore 
totally consistent with observations. 

III. VACUUM ENERGY AND ULTRA-LIGHT 
SCALARS 

The quadratic divergence proportional to i7^(t)A^ has 
the interesting property that, if we set the cutoff Ag at a 
value of the order of the Planck mass Mpi, it becomes of 
the order of the critical density of the Universe at time 
t, since pc{t) = (3/87r)iJ^(t)Afp[. This is in sharp con- 
trast with the quartic divergence A^i, that for Ac ^ Mp\ is 
more than 120 order of magnitudes larger than the obser- 
vationally allowed value. Thus, the quadratic divergence 
might offer a hint for understanding the coincidence prob- 
lem, namely why the vacuum energy density that drives 
the acceleration of the Universe is of the order of the 
critical density today. In this naive form, however, this 
argument does not lead us very far, since the quadratic 
divergence is simply reabsorbed into a renormalization of 
Newton's constant and docs not leave observable effects. 
Still, the fact that a term proportional to H'^{t) in the 
vacuum energy density might naturally be of the order 
of the critical density is quite intriguing and stimulates 
further investigation. In this section we will see that it 
is indeed possible to prevent the term cx H'^{t) of being 
reabsorbed into a renormalization of Newton's constant. 
However, to do so we are required to introduce an ultra- 
light degree of freedom into the theory. In this paper, 
for definiteness, we will consider a scalar field. However, 
most of our considerations in this section go through as 
well if we introduce an ultra-light massive graviton in- 
stead. The role of the scalar field is then played by the 



9 



helicity-0 component of the massive graviton. Ultra-light 
scalars with m ^ Hq where first discussed, in the context 
of pseudo-Nambu-Goldstone bosons, in [32]. In ref . [33] it 
was first suggested that quantum effects from ultra-light 
scalar fields could lead to late-time acceleration. 



A. Ultra-light scalars and the efFective action 

We begin by investigating the most general structure of 
the renormalized VEV of the energy-momentum tensor. 
First of all, to avoid misunderstandings, let us empha- 
size that the answer to this question does not depend on 
the regularization scheme employed. This may affect the 
structure of the bare (0]r^,y]0), but not that of the renor- 
malized one. In a theory without ultra-light scalar fields, 
a generally covariant regularization gives the most gen- 
eral structure of the bare (0]T^y]0) as in eqs. (1) and (10). 
In this case even the counter-terms are generally covari- 
ant, so the renormalized (0)7^,^10) can in general have a 
term proportional to g^i^ from the variation of the volume 
term (which is set to zero if one adopts the proposal in 
section II C, or otherwise reabsorbed into a renormaliza- 
tion of the cosmological constant). Further it can have 
a term proportional to the Einstein tensor G^i, from the 
variation of the Einstein-Hilbert term, reabsorbed into 
Newton's constant, and contributions from the variation 
of the terms quadratic in the Riemann tensor. On di- 
mensional grounds the latter are suppressed by a factor 
(E/Mpi)"^ with respect to the Einstein-Hilbert term at 
the energy (or curvature) scale E, and can therefore only 
be relevant at Planckian energies. 

As we have discussed, using a momentum-space cut- 
off slightly complicates the analysis as the explicit cal- 
culation of the bare (0|Tjuy|0) in such a scheme produces 
non-covariant contributions. These are compensated by 
introducing non-covariant counter-terms, chosen so to re- 
store general covariance for the renormalized quantities. 
Thus, in the end the renormalized (OjT^^lO) is indepen- 
dent of the regularization scheme, as it should. In the 
following, we will use a generally covariant regulariza- 
tion, so the general structure of the bare (0]T]ui^]0) is the 
same as that of the renormalized (OjT^i/jO). 

To understand if it is possible to obtain a more general 
result for (OjT^^lO) it is useful to recall that the VEV of 
T^i, can be obtained by taking the functional derivative 
of the effective action with respect to the metric, see e.g. 
[14, 15]. The effective action Sca[ginA is derived from the 
matter action 5^ [ff/i^ , f/'] by integrating over the matter 
fields, here collectively denoted by tjj 



iSeff [g/j, 



-I 



(33) 



The matter energy-momentum tensor is defined as usual 
by T^^ = -(2/7=5) SSm/Sg^"". In the path-integral for- 
mulation its vacuum expectation value is given by 



(0|7;.|0) = 



iSmigp^^M 



(34) 



and can be rewritten as 



-2 



2 SS,e 



-9 Sgt'^ 



(35) 



The effective action defined in eq. (33) is the relevant 
quantity for describing the physics at an energy scale 
that is sufficiently low compared to the mass of all mat- 
ter fields, so that the matter fields can be integrated out. 
In FRW the relevant scale is given by the Hubble param- 
eter H{t) at the time t of interest. Thus, one is implicitly 
assuming that all matter fields are very massive with re- 
spect to the Hubble scale, m ^ H(t). However, if the 
theory includes a scalar field with mass m<H{t), this 
assumption is violated since such a field cannot be inte- 
grated out completely. In particular, in the presence of 
an ultra-light scalar with m ~ /fo, the assumption is vi- 
olated for all times t < to (where to denotes the present 
age of the Universe). More precisely, at time t we can 
integrate out the Fourier modes (pk (where k denotes the 
comoving momentum, related to the physical momentum 
by kphys = k/a(t)) whose frequency satisfies 



Wfe = + (A;Va2) > H{t) . 



(36) 



The low- frequency modes with k such that ujk ^ H{t) re- 
main in the effective action and should be treated as a 
classical background, on the same footing as the FRW 
metric .9^,/.** It is then convenient to split the field 
into its low- and high-frequency parts. We first expand 
^(x) in Fourier modes, 



<t>{x) 



I (27r 



(27r)3V2fc 



and define 



(j){x) = x{x) + (p{x) 



(37) 



(38) 



where in x(x) the integral over (Pk runs only over the 
comoving momenta with ]k] < A;«(t) while in ip{x) it 
runs only over comoving momenta with |k| > k^,{t). The 
value k^ that separates the two regimes can be defined 
by rewriting eq. (36) as ujk, = vH{t) where is a dimen- 
sionless constant with v^l.^ Therefore 



k^{t) = a{t)^/v^m{t)-m? . 



(39) 



® Our approach here differs from that in [34] , where the authors in- 
troduce an ultralight scalar field with m ~ -ffoi and then derive 
an effective action for the gravitational field only, by integrat- 
ing out all modes of the scalar field, including its low-frequency 

modes. 

® As usual in this separation the constant u is in principle arbitrary, 
and the evolution of the low-energy effective theory with respect 
to V is described by renormalization group equations. 
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Figure 1: A Feyman diagram that contributes to the gener- 
ation of a term f{x)R in the effective action. The wavy hne 
represents the classical gravitational field, the solid line in the 

loop is the field ip and the dashed line represents the field x- 
A graph with n external x lines contributes a term oc x"^^ to 
the effective action (here n = 4). 



Of course, the low-momentum modes exist only if m and 
t are such that m < uH{t). At times where m ^ H{t), 
all modes become heavy and cj) can be integrated out 
completely, or equivalcntly fc^ (t) = 0. Observe that the 
separation between x and is time-dependent. Denoting 
all other massive matter fields collectively by if), the effec- 
tive action now becomes a functional of both the metric 
and X) given by 



/ 



(40) 



In principle S^filQuv, x] will include all possible terms that 

are consistent with general covariance. 

In particular, there can be a term proportional to 
J d^XyJ—g /(x)-R; where /(x) is some function whose ex- 
plicit form depends on the details of the original action. 
Diagrammatically, this term comes from graphs such as 
that in fig. 1, where the wavy lines represent the back- 
ground gravitational field, the solid line circulating in 
the loop represents the high-frequency modes and the 
external dashed lines represent x- A graph with the in- 
sertion of n external x fines contributes a term propor- 
tional to x^-R to the effective action. Observe that the 
X field only occurs in external lines, since it is not inte- 
grated over. Further notice that a graph with two exter- 
nal wavy lines produces the term in i? quadratic in the 
background metric. Terms cubic and of higher order in 
the background metric coming from R are obtained from 
graphs with three or more external wavy lines. The inte- 
gral in the loop runs only over high-frequency momenta, 
i.e. fc, < < 00, and can be regularized in the UV 
using e.g. dimensional regularization. 

The graphs shown in fig. 1 assumes the existence of a 
interaction term in the potential 'V{4>)- After splitting 
(/) = X + </' this produces a vertex X'W- With a A(/)^ term 
in the potential we had vertices such as XXW instead. 



where two external x fines arc attached at each vertex 
in the loop. Furthermore, the gravitational field must be 
similarly split as — g^iv-^^ixv, where the background 
g^iv contains the low frequency modes, and basically co- 
incide with the background FRW metric, while ft^j/ con- 
tains the high-frequc!ncy modes and describes quantum 
fluctuations over the background. In eq. (40) we must 
therefore integrate also over /i^^,, so there are also graphs 
such as that in fig. 1 in which the internal (p loop is re- 
placed by a loop where hp,v circulates. The non-linear 
terms in the gravitational action provide vertices such as 
ghh that allow us to attach external gy,y lines to the /i^j^ 
loop. Similarly the interaction of </> with gravity provides 
vertices such as x^^ that allows us to attach external 
X fines to the loop. Therefore terms of the form f{x)R 
are generated through the gravitational interaction even 
if V{(j)) = (l/2)m^(/)^ and there is no direct interaction 
between x and (p. All these graphs are generically non- 
zero, and therefore a term f{x)R is in general generated 
in the effective action, even if there was no term f{(j))R 
in the original action for </>. 

Among all terms that can be generated by the integra- 
tion over the high-frequency modes, at the phenomeno- 
logical level the term f{x)R is the most interesting one, 
since it involves the lowest number of derivatives. There- 
fore it is not suppressed at scales below the Planck mass, 
contrary to terms involving f?^, R^i.Rf^", etc. or interac- 
tion terms involving derivatives of x- Therefore the most 
general effective action, for generating a cosmologically 
relevant (0|T^iy|0), is of the form 




(41) 



167rG 



-yd^Xd.X-V{x) 



plus higher-derivative operators, that are not relevant at 
scale much below the Planck scale, and possibly non- 
minimal kinetic terms, as in k-essence models. 

The integration over the high-frequency modes in FRW 
has also been studied, in a somewhat different language, 
in the context of stochastic infiation [35-37]. For in- 
stance, for a minimally coupled scalar field ^ in De Sit- 
ter space, the result for Ses[gi^v,x] can be written in the 
form [36] 

where is an action from which one derives the classical 
equation of motion for Xi in the form 

X + SHx- a-^V\ + V'{x)=i + ■ (43) 

Here ^ is a random field with a two-point correlator 

{ax)ax')) = A{x - x') , (44) 

and P[^] in eq. (42) is given by 

P[^]=exp|-^ J d'^xd'^x' ^{x)A-\x,x')^{x')^ . (45) 
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Defining again the high-frequency modes as the modes 
with uJk > vH ^ the form of A{x — x') in the most in- 
teresting limit v>Ji has not been computed, while in the 
limit 1/ -C 1 it is given by [36] 

A{x -x') = -^ '- S{t - t') + 0{v^) , 46 

47r^ uH ar 

where r = |x — x'|. The above formulation is physically 
quite intuitive as it allows us to interpret X as a classical 
field driven by some stochastic noise. This is the typical 
picture of stochastic inflation. If one was to integrate 
out ^ from eq. (42) one would find the effective action 
Ses[gij,v,x] that can be directly obtained from diagrams 
such as that in fig. 1. 

It is also instructive to directly compute (0|rp,y|0) 
based on the sum of zero-point energies, and to see how 
a term f{x)R could be generated in the effective action 
even if it was not present at first. In the usual setting 
with only massive fields, e.g. a massive scalar field $, one 
expands the field in terms of creation and annihilation 
operators 



(27r)V2A; 



ak$fc(i)e 



ik-x 



-?k-x 



(47) 

The modes (t) are obtained by solving the equation of 
motion in the given FRW metric, here 



< + 2-$'fe + fc2$, = o, 



(48) 



where the prime denotes derivatives with respect to con- 
formal time T]. Then the vacuum energy density is given 

by 



1 f (Pk 
=2/(2^ 



(49) 



where some regularization of the momentum integral is 
understood, e.g. a cutoff over momenta or dimensional 
regularization. Thus, according to the standard rules of 
semiclassical gravity, the field $ is treated as a quantum 
field, while the metric g^j^ is treated classically, and the 
modes of the quantum field are determined solving an 
equation in this classical background. The energy den- 
sity (0|Too|0) calculated in this way depends explicitly on 
the Hubble parameter H(t) because the modes $fe(?7) de- 
pend on the background metric, in this case on the scale 
factor a (77) and conformal time \r}\ ^ {aH)~^ of FRW. 
For example, in a De Sitter epoch or during matter dom- 
inance (MD), 



1 



a(r?) 



kri 



-ikrj 



(50) 



In the limit fcry — > 00, and therefore fcphys = k/a ^ H, 
these modes reduce to the usual plane waves of flat space 
(this corresponds to the most natural choice of vacuum 
state, namely the Bunch-Davies vacuum; the most gen- 
eral form of the modes is a superposition of positive 



and negative frequency modes with arbitrary Bogoliubov 
coefficients). Therefore the leading UV divergence in 
eq. (49) is the same as in flat space, i.e. the term if we 
regularize with a cutoff in momentum space. However, 
the term i/kr] in eq. (50), as well as the time dependence 
in the factor 1/0(77), affect the subleading divergence, and 
produce a quadratic divergence oc H'^{t)K^. (For details 
on these calculations see e.g. ref. [10].) 

When there is also an ultra-light scalar field (f) the 
above scenario changes, because </> must be split into its 
low- frequency modes x and its high-frequency modes 93. 
The field x should be treated as classical on the same 
footing as the metric g^j^i, since, on physical grounds, 
modes with wavelength larger than the horizon are ex- 
pected to undergo decoherence. (See ref. [37] for a more 
careful discussion of this point.) The high-frequency 
modes should instead be quantized in the classical 
background, given both by and x- In particular, 
the equation of motion for Lpk that replaces eq. (48) will 
depend both on the scale factor and on the homogeneous 
field x(^)- It is derived from the equation of motion of 
(f) by setting (f) = x + <f and expanding ip in momen- 
tum modes (pk- As a result the modes (50) will acquire 
a dependence on x which will generically translate into 
a x-dependence of (0|T^j,|0). In general, this produces 
all terms which are consistent with the symmetries of 
the problem, and therefore also terms corresponding to 
f{x)R in the effective action. 

Finally, we observe that similar considerations hold 
for massless conformally- coupled fields, including the 
electromagnetic field. If it were not for the conformal 
anomaly, in FRW (which is conformally equivalent to 
Minkowski space) these fields would decouple from the 
metric and could be integrated out trivially. Because of 
the conformal anomaly at the one-loop level these fields 
do not decouple but still they can be integrated out ex- 
actly, giving rise to the anomaly-induced effective action. 
However the fact that, technically, it is possible to inte- 
grate them out exactly does not mean that, physically, 
this is the correct thing to do. In any case, the long- 
wavelength modes of these fields should be kept in the 
effective action, on the same footing as the classical met- 
ric, and produce a classical background (whether of a 
scalar field, of an electromagnetic field, etc.) which in 
general will be relevant on cosmological scales. Excluding 
from the integration the modes with physical momenta 
fcphys < kt,/a{t) ~ H{t) could also in principle affect the 
non-local terms of the anomaly-induced effective action. 



B. Interacting vacuum fluctuations 

We next discuss the properties of (0|Tpi,|0) derived 
from the effective action (41). Models of this form have 
been studied in the literature as generalizations of the 
simplest quintessence models [38-43]. Of course such 
models can be introduced and investigated in their own 
right, without being related to zero-point fluctuations. 
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Our line of arguments leading to such models, based on 
the most general form of (0|r^^|0) in the presence of 
ultra-light scalar fields, gives further motivation for their 
study. This insight will also allow us the recast the cos- 
mological equations of the model in an interesting form, 
which can be seen to describe the interaction between 
vacuum fluctuation of the field 4> and the classical part x 
of the same field. 

Specializing to a FRW metric and to a spatially homo- 
geneous field X = x{t)i tbe equations of motion derived 
from the action (41) are [38-43] 



2H + 3iJ2 
X + 3Hx + V'{x) 



SttG 
3 

~8nG 



{PB + Px), (51) 
Px) > (52) 
{H + 2H')f'{x), (53) 



where Pb, Pb are the energy density and pressure of the 

matter-radiation background, respectively, and p^, 
are those associated to x as given below. The prime 
denotes the derivative with respect to x- Denoting the 
total EOS parameter by t«tot = {Pb +Px)/ipB + Px)' 
two Priedmann equations are combined to find 



H = --H^{l + Wtot), 

such that eq. (53) can be recast as 

3 



(54) 



X + 3Hx + V'{x) 



H^{l-3wtot)f{x)- (55) 



It is convenient to separate py^ and p^^ into two parts, 

Px= Px+Pz, Px=Px+Pz, (56) 

where 



Px = ^x' + V{x) , 



1 .2 

Px=^X 



V{x) , (57) 



are the standard contributions of the classical part, while 



Pz = 



8nG 
1 

Pz = - 



H^f{x) + Hf{x) 



(58) 



SttG 



i2H + 3H')f{x) + 2Hf{x) + fix) -(59) 



Using again eq. (54), eq. (59) can also be rewritten as 
1 



Pz 



8nG 



3wtotH'f{x) - 2Hf{x) - fix) 



(60) 



The subscript Z, which stands for "zero-point fluctua- 
tions" , stresses that in our approach pz and pz can be 
interpreted as a contribution coming from vacuum fluc- 
tuations, i.e. from the VEV (0|T^^|0) obtained through 
eq. (1) from an effective action that includes the contri- 
bution f{x)R which generically appears in the presence 
of the ultra-light field %• 

A rather interesting property of pz is that it is of the 
order of the critical density at all times. For instance. 



recalling that the critical density at time t is pdt) = 
i3/8TrG)H^it), we sec that the term i3/8TrG)H^{t)fix) 
in pz is just /(x)Pc(O- 1^ the slow-roll regime x(t) is 
roughly constant, and this term provides an approxi- 
mately fixed fraction of the critical energy density. In 
this way we obtain a term of order Mpj7J^(t) in the en- 
ergy density, which is not reabsorbed into a renormaliza- 
tion of Newton's constant because its prefactor contains 
the function fix)- 

The above equations can be rewritten in a form which 
allows for an appealing physical interpretation, namely 
that the "classical" energy-momentum tensor described 
by Px and px interacts with the vacuum fluctuations, 
described by pz and pz- The interaction is governed 
by a coupled energy-momentum conservation equation 
(thereby providing an explicit realization of a mechanism 
proposed in [10, 11]). To this purpose, it is convenient 
to transform the time derivatives into derivatives with 
respect to a; = In a, which we denote by d^. By means of 
eq. (54) we can then write 



pz = Pcif + dxf) , 
Pz = wzpz , 

(1 + 3wtot)d^f + 2dlf 



Wz = Wtot - 



Qif + d.f) 



(61) 
(62) 

(63) 



If the field x does not exchange energy with the other 
forms of matter included in the background action Sb (as 

we have indeed already assumed when deriving eq. (53)), 
its energy density and pressure satisfy the conservation 
equation 



Px + -iHiPx+Px) =0. 
Splitting p^ and p^^ as in eq. (56) and defining 

_px ^ (l/2)x^ ~ Vix) 
- PX (l/2)x2 + Vix) ' 

we can write the continuity equation as 

Pz +3Hil + wz)pz = Q, 
Px +3Hil + ivx)px = -Q, 



(64) 



(65) 



(66) 
(67) 



where Q(t) is the energy density transfer rate from X to 
Z. Inserting the explicit expressions for pz and wz given 
in eqs. (61) and (63) we find 



Q 



-il-3w^,t)Hit)p,it)dJ. 



(68) 



Since pdt) and Hit) provide the natural scales for the 
energy density and for the inverse of time, respectively, a 
natural scale for the rate of energy transfer is Hit) pdt). 
We sec from eq. (68) that, with respect to this scale, the 
actual rate Q differs by a factor proportional to dxf and 
to (1 — 3wtot)- Since Hd^f = /, we see that the actual 
rate of energy transfer is governed by the time variation of 
f{x), and therefore in a slow-roll phase \Q\ <C H{t)pdt)- 
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Furthermore, during RD i^tot ~ providing a further 
suppression factor. However, note that Q does in general 
not vanish exactly in RD because of the contributions 
from pz and matter. We will see this for a typical exam- 
ple in the next section. 

Conservation equations with energy exchange such as 
eqs. (66) and (67), that even go back to Bronstein (1933) 
(see [5]), are at the basis of much recent work on inter- 
acting dark energy/matter models (sec e.g. [39, 44-49] 
and references therein). Typically the functional form of 
Q{t) is chosen on purely phenomenological grounds. It is 
interesting to see that in this model the form of Q{t) is 
explicitly predicted in terms of dxf -^^ 

In summary, the energy-momentum tensor of the clas- 
sical part, r^j^ = [pxi o?^ijPx) a-iid that of the vacuum 



fluctuations (0|r^i.|0) 
separately conserved: 



{pz, a 6^jpz) are in general not 



V'*(0|T^.|0) = -V''T0,^0. 



(69) 



Rather, they interact through a coupled energy- 
momentum conservation equation. 



C. Cosmological evolution and tracking dark 
energy 

The cosmological evolution obtained from eqs. (51)- 
(53) has been studied in detail, see e.g. refs. [38-43]. To 
illustrate the typical evolution, we consider a model with 



fix) 

16ttG 



Fix) = ^Cx' 



2 2 



(70) 



and we take for definiteness ^ > 0. It is useful to intro- 
duce the dimensionless parameter 



(71) 



The ratio of the term 6H^F in pz to the potential V 
is equal to / {^j? Hq) . Therefore, if /z = 3, these two 
terms become comparable at the present epoch, while in 
the past 6H^F ^ V. For larger values of p the potential 
V begins to give a significant contribution to the energy 
density somewhat earlier than today, while for smaller p. 
it will only become important in the future. The evolu- 
tion of X, determined by eq. (55), is different depending 
on whether V'{x) or 3H^ {l—3wtot)F' ix) dominates. The 
ratio of these two terms is 



3^2(1 _ 3«;tot)i^'(x) 3if2(i _ swtot) 



V'ix) 



p?Hl 



(72) 



To make the connection with the discussion just after eq. (15) it 
is useful to introduce the EOS parameter defined by the 

condition pz + 3H{1 + Wt)pz = 0. Using eqs. (61) and (54) we 
get = «>tot - (d:cf + d'^f)/[3(f + d:cf)]- Observe that d^f = 
implies Wt = lutot consistent with eq. (15). Comparing with 
eq. (63) wc sec that tuz = ui, +tD where m) = [l—3wtot]dxf/[6{f+ 
dxf)] and therefore in general u>z Wa,. 



Since obscrvationally wj^^j w —0.7 today, this ratio be- 
comes of order unity at the present epoch for p, ~ 3. In 
the past it was much larger than one. This is true even 
in RD despite the fact that, during RD, (1 — 3wtot) 1- 
To see this in a simpler setting, lets us at first in- 
clude just the contribution of matter to Wtot in the RD 
phase. Including both matter and radiation, deep in RD 
where pM *C Pr and the scale factor a <^ a^q (where 
aeq = fl%/fl°M ~ 3 X 10"^) we have 



M'tot 



(l/3)Pfl 1 1 



PR + PM 3 1 -h a/fle 



Therefore 



(1 - 3w,ot)H^ ^ —H^ . 



(73) 



(74) 



In RD a cx tM^ and H"^ (xl/t^, so (1 - 3wtoi)H^ fx t-^/^ 
and, going backward in time, the growth of more than 
compensate the decrease of (1 — 3ti;tot)- 

More importantly, in this model there is a significant 
contribution from DE to Wtot even in RD. This contri- 
bution, contrary to that of matter, is not parametrically 
suppressed in the limit a/ acq ^ 1. To compute it in the 
matter and radiation era we can safely neglect px ■ Then 
we find from the expressions in the last section 



1 - Swt. 



(1 - / - d,f)/{l + aeq/g) + d^f + dlf 
1 - / - dJ/2 



(75) 

which correctly reduces to eq. (73) when / = const. 
Thus, the combination (1 — 3w7iot) is gcnerically non- 
zero, and for p = 0(1) the evolution of x during RD and 
most of MD is dominated by the non-minimal coupling 
F{x)R- The potential V{x) only becomes relevant when 
approaching the current epoch, while for larger values of 
p the potential becomes relevant earlier. Eventually, the 
field X will exit the slow-roll phase and perform damped 
oscillations around the minimum of V{x) at % = 0. 

In this model the total DE density px + pz tracks 
the total energy density pdt), as we can see from fig. 2. 
This is basically due to the combination of two factors: 
dtiring RD and deep into MD the tracking is ensured 
by the term /(x)Pc(i) in Pz- This term, however, can- 
not be the dominant one today, since it is proportional 
to H'^it), and in the recent cosmological epoch a time 
dependence poE H^it) is obscrvationally excluded. 
Thus, approaching the recent epoch, we must switch to a 
regime where px dominates and x is still approximately 
slowly-rolling such that Wx is close to —1, see eq. (65), 
and the DE density is approximately constant. 

With suitable values of the parameters ^, p and of the 
initial conditions we can get values for Ode and wde 
today of the order of the obscrvationally preferred ones. 
For instance, in fig. 2 we used p = 24.3, ^ = 5 x 10"'' 
and we set the initial conditions Xin = 4 and xL = 
deep in RD (measuring x in units SttG = 1). With these 
values we get Ode — 0.73 and wde = —0.952 today. 
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Figure 2: Left: logp for radiation (red, dotted), matter (blue, dot-dashed) and dark energy (black solid) as functions of log a. 
Right: the same plot on a linear vertical scale, in terms of fli — pi{t) / pdt), with i = R, M, DE. We used the values pb = 24.3, 
^ = 5 X 10~* and the initial conditions Xin ~ 4 (in units 8itG — 1), Xin = deep in RD. 




Figure 3: Left: Wtot (red, solid) and wde (blue, dashed) as functions of log a for the same parameters as in fig. 2. Right: the 
separate contributions fix = px(t) / Ptotif) (black, dashed) and Q.z ~ pzif) / ptotif) (brown, dotted). 



The evolution of the EOS parameter wde = Pde/pde 
is shown in the left panel of fig. 3, together with the 
evolution of Wtot ■ This shows that a realistic cosmological 
evolution can be obtained in this model. 

It should however be stressed that the coincidence 
problem is not really solved by the fact that the DE 
density tracks the total energy. The total dark energy 
density pde = Px + Pz switches from an early regime 
where it is dominated by pz, and hence evolves approxi- 
mately as H^{t), to a late regime where it is dominated 
by Px, and therefore is approximately constant. The 
transition between these regimes can nicely be observed 
when looking at the separate evolution of pxit) / pdt) 
and pz{t)/ Pc{t), plotted in the right panel of fig. 3 for 



Furthermore, in this model the limits on the time variation of 
Newton's constant require 32ttG{^xo)^ < ^bO' where xo is the 
present value of x{t) and ujbd is the Brans-Dicke parameter [41]. 
Using the bound ujbd > 4.2 X 10^ inferred by the Cassini mission 
[50] and measuring x in units SttG = 1, this bound reads ^xo < 
2.4 X 10~^. This is satisfied by the solution shown in fig. 2, for 
which 5x0 ^ 1-9 X 10-3. 



the same model parameters used for the other figures. 
The time when the transition takes place is controlled by 
ratio between the mass of x and the Hubble constant. 
In fact, it is the parameter p, = m/ {S}/"^ Hq) that is cru- 
cial here: the larger p the earlier the transition. As we 
argued above, ii p — i the transition happens about at 
the current epoch. Notice though, compared to ACDIVI, 
the coincidence problem is still sensibly alleviated. This 
is clearly seen from fig. 4, where the ratio PDE(i)/ptot(i), 
computed in the model with F{x) and V{x) given in 
eq. (70), is compared to the same ratio in ACDM. 



Exploring the parameter space to a certain extent, we 
find that the figures shown represent the typical behavior 
of this model, as long reasonable initial conditions are 
set. In principle, a deeper study of the stability and the 
parameter space of the system would be needed, but it is 
not the aim of this work to do so. The reason is the lack 
of motivation for the specific choices for F{x) and V{x) 
at this point. Therefore we will discuss a more general if 
slightly approximate approach in the following. 
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Figure 4: The ratio Q,-DE{t) = poKit) / ptotit) , computed in 
the model with F(x) a^nd V{x) given in eq. (70) (solid, black) 
compared to the ratio O.A{i) = ph/ ptot{t) in ACDM (dashed, 
red). 



D. Comparison with observations 

The model that we have discussed has a DE den- 
sity pDE(i) = Pxit) + pz{t), where px{t) is the typi- 
cal quintessence energy density due to the rolling of the 
scalar field x in the potential V{x), while pz[t) comes 
from the term F{x)R in the effective action. With a suit- 
able choice of the mass m in V{x) the term px comes to 
dominate near the present epoch, while it is totally neg- 
ligible at early times, and its EOS parameter wx is close 
to —1 today, as required by the observations. In con- 
trast, pz{t) is proportional to H^{t). In order to have a 
viable model, pz must be subleading today since a DE 
term proportional to H^{t) is observationally excluded 
(see e.g. [47-49]). However, the fact that this term is 
approximately a constant fraction of pdt) means that its 
integrated effect over the whole history of the Universe 
could give observable effects. 



1. A more general parametrization of pj^,^ 

The aim of this section is to perform a comparison of 
such a model with cosmological observations, expanding 
on the results that we have presented in [11]. A possi- 
ble way to proceed could be to choose a particular form 
of the functions F{x) and T^(x), such as those given in 
eq. (70), and find the corresponding limits on the param- 
eters ^ and m. However, such limits would be tied to the 
specific choice made for the functions F{x) and V{x)- In 
the absence of strong theoretical motivations for any such 
choice, it seems more useful to extract some general fea- 
tures of the DE sector of such models, and parametrize 
the DE density in a form that might have a more general 
validity. This can be done by recasting eqs. (57) and (58) 
in the form 



PDE = + 



\ (dxxf 



6iF + d,F) 



(76) 



where we have again transformed the derivative with re- 
spect to t into the derivative with respect to a; = In a, 
denoted by d^- By means of the Friedmann equation, 
H"^ — {8nG/3){pii + pm + pde), and the definition 



8ttG 



\ {d^xf + 6(F -f d.,F) 



we can write the DE density in the form 

V + {pB, + Pm)Z 



1 



(77) 



(78) 



Finally, we implicitly define the three functions WQ{a)^ 
eii{a) and Em (a) and the constant C by writing 

PDE - Ca'^^^+^»^'''>'>+pR{a)eR{a)+pM{a)eM{a), (79) 
such that 

Z 



ei?(a) = eA/(a) = e(a) 



wo[a) 



-1 - ^ In 
3 



V 

cz 



1 



e{a) 



(80) 
(81) 



When Wo is constant, the term Ca^'^'^^^^"^ in eq. (79) 
corresponds to the usual dark energy density of a minimal 
quintessence model. On the other hand, the presence of 
the terms pReu and pm^m clearly displays the tracking 
behavior of the DE density. Choosing a specific model 
for the functions F{x) and V{x) corresponds to making 
a choice for wo{a) and e(a). The constant C is fixed by 
the DE density today, or better, by a combination of the 
model parameters in F and V and the initial value of x- 
Going back to the example discussed in the last sec- 
tion, F{x) = (l/2)^x^ and V{x) = {l/2)m'^x'^ , we solve 
eqs. (51) and (53) numerically and show the evolution of 
the functions wo{a) and e(a) in fig. 5. We see that in the 
recent epoch where the term Ca~'^'-^+"'°^ is relevant wa 
can roughly be approximated by a constant, within a few 
percent. Moreover, the function e{a) is constant to excel- 
lent approximation in RD and decreases linearly with In a 
in MD. This is the reason why we separated the terms 
into prEr and pm^m, as either one will be negligible 
when the other one is dominating. Although e(a) raises 
again when the present epoch is approached, see fig. 5, at 
this stage the contribution from pm^m in the DE density 
is subdominant compared to the term Ca~^^^+'^'') com- 
ing from Px (see the right panel in fig. 3). Accordingly, 
the behavior of e{x) in this regime is not very relevant. 
Therefore, we can approximately write 



Em (a) 



ej\,/ + _B In a 



(82) 
(83) 



with independent constant parameters and cm- The 
constant B, on the other hand, is fixed by requiring 
the two branches to coincide at equality: B = {er — 
Em)/ In acq. 
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Figure 5: The evolution of the functions e = — Z] and wq defined in eq. (79) are shown for the model and parameters of 
Sec. IIIC, see fig. 2. In the left panel we also show the en — em approximation with en — 0.00806 and ea/ = 0.00793 (blue, 
dashed) . 



The fact that can be taken as constant is a conse- 
quence of X rolling very slowly during RD, in fact it does 
essentially not move at all, see also ref. [43]. Thus, the 
function Z{x) is virtually constant. The fact that in the 
MD era e evolves linearly with In a is because x begins 
to roll, still slowly though, and therefore with In a as this 
is the natural variable for cosmological evolution. Thus, 
this behavior appears to be quite general for typical vi- 
able models, and is likely to be not restricted to the spe- 
cific model F{x) = (l/2)^x^ and V{x) = il/2)m'^x'^ that 
we have studied here. We will therefore use eq. (79) as 
our basic expression for the DE density, with u'o(a) = wq 
and e(a) given in terms of the constans en and sm- 
This model then is a generalization of the one studied 
in rcf. [11], which corresponds to the limit e/j = ejv/- It 
is a two-parameter generalization of wCDM, which in 
turn corresponds to the limit e/? = eM = 0. We will 
generically use the label wZCDM for this generalization 
(including the limiting case tR = studied in [11]). 

2. Results 

We estimate the model parameters with a full like- 
lihood analysis using modified versions of the publicly 
available Markov-Chain Monte-Carlo sampling code Cos- 
moMC [51, 52] and the CMB-Boltzmann code CAMB 
[53, 54]. We assume a flat FRW background and put a 
prior on the age of the Universe to be between 10 and 20 
Gyrs. Then we use the following datasets: constraints on 
the current Hubble parameter from Hubble Space Tele- 
scope (HST) observations [55] and on the number of rel- 
ativistic degrees of freedom at Big-Bang Nucleosynthesis 
(BBN) [56]; the Union2 Compilation of Type la super- 
novae (SNe) of the Supernova Cosmology Project [57]; 
the angular power spectra of temperature and polariza- 
tion anisotropics in the cosmic microwave background 
(CMB) from WMAP7 [58], ACBAR [59] and CBI [60]; 
and the Baryon Acoustic Oscillations (BAO) data from 
the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) 



[61]. 

We perform the analysis for the four following models: 
in wZCDM2 we vary {lUo, e^j, ea/}; in wZCDMl we vary 
{wq, e} and set e/j = cm = e; in luCDM we fix e = in 
addition and, finally, in ACDM we fix = — 1 as well. 
On top of these DE parameters we vary the common set 
of six base parameters. These are the slope Ug and am- 
plitude log[10^°As] of the spectrum of primordial scalar 
curvature perturbations (modeled as a power law nor- 
malized at fc = 0.05 Mpc~^), the depth to re-ionization 
T, the physical baryon energy fraction flbh^ , the physical 
CDM energy fraction ^Ij^uh^, and 100 times the ratio of 
the sound horizon to the angular diameter distance to 
the last-scattering surface, 9. We use flat priors for all 
parameters, set adiabatic initial conditions for the evolu- 
tion of the cosmological perturbations and ignore vector 
and tensor modes for simplicity. 

As opposed to the cosmological constant A, in the case 
of dynamical DE models we need to model the evolution 
of its inhomogeneities. If we were interested in a specific 
model, we should solve the perturbation equations de- 
rived from its Lagrangian. However, as discussed in Sec- 
tion HID 1, in the context of this work we are not really 
interested in analyzing a model with a specific form of 
the functions F{x) and V{x)- We rather wish to extract 
more general features of a class of DE models suggested 
by our theoretical considerations. In this spirit, we will 
at first assume the DE to be a perfect fluid with a con- 
stant rest frame sound speed Cs = 1, as it is usually done 
in the wCDM model, and we will later check how our 
results depend on this assumption. 

The results are summarized in table I. We show 
the marginalized posterior probability distributions for 
wZCDMl and it;ZCDM2 in fig. 6, comparing both cases 
with wCDM and ACDM. For wZCDMl we find lower 
and upper limits of 

-1.25 < Wo < -0.908 , -0.0201 < e < 0.0460 , 

(84) 

at 95% C.L. For wZCDM2 we find lower and upper limits 
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ACDM 


wCDM 


wiZCDMl 


w;ZCDM2 


(...) 


(...) 


(...) 


(...) 





0.966 


0.011 


0.963 


0.013 


0.967 


0.014 


0.971 


0.016 


log[10i°^.] 


3.084 


0.032 


3.088 


0.032 


3.087 


0.032 


3.084 


0.034 


T 


0.087 


0.014 


0.086 


0.014 


0.087 


0.015 


0.089 


0.015 




0.023 


0.001 


0.023 


0.001 


0.0226 


0.001 


0.0225 


0.001 




0.113 


0.003 


0.115 


0.004 


0.115 


0.005 


0.111 


0.007 


e 


1.040 


0.002 


1.040 


0.002 


1.039 


0.003 


1.042 


0.006 


Wo 






-1.043 


0.079 


-1.067 


0.087 


-1.11 


0.11 


€ 










0.0106 


0.017 




















0.041 


0.048 
















-0.0068 


0.026 



Table I: Means and standard deviations of the marginalized likelihoods. We observe that all additional DE parameters are 
consistent with ACDM at the 0.68% C.L. for all models wCDM, wZCDMl as well as wZCDM2. All models have unit rest 
frame sound speed, Cs = 1. 



of 

-1.35 < Wo < -0.903 , -0.0465 < cm < 0.140 , 
-0.0558 <€R< 0.0446 , (85) 

again at 95% C.L. In general, we observe a slight pref- 
erence for models with Wq < —1 and Cm > 0, i.e. early 
DE models with a recent phantom-crossing. However, all 
new parameters are consistent with their ACDM values 
at the 95% C.L. as this is also the case for ?i,'o hi w;CDM. 
The results are consistent with those for other early DE 
models found in the literature, e.g. [62-64]. However, 
note that our model has a slightly different evolution 
during radiation domination than the widely used pa- 
rameterization proposed in ref. [65], and usually for that 
model phantom-crossing is excluded with the hard prior 
wq> —1 which we do not employ here. 

In fig. 7 we show the combined likelihood contours for 
a subset of the parameters of wZCDMl and compare 
the constraints with those for the wCDM model. The 
slight degeneracy in the (e, wq) plane comes from the 
fact that ■u;o(a) depends logarithmically on e(a) as can be 
seen from eq. (80). For positive e the data prefers wq < 
— 1. Secondly, we note that the constraints on Ode are 
virtually unchanged by introducing the new parameter e, 
as there is basically no degeneracy in the (e, ^de) plane 
and the marginalized posterior of for J^de is unchanged 
as well. Finally, e is mildly correlated with the spectral 
index, Ug, as our analysis shows. This leads to a slight 
shift of the posterior to larger values of rig in wZCDMl 
as compared to wCDM. 

In fig. 8 we show the combined likelihood contours 
for wZCDM2, comparing to wCDM. The correlation be- 
tween Cm and Wq persists, which means that e in the 
simplified wZCDMl model is more or less reflecting eM- 
This is because the observations of the SNe and BAO are 
more sensitive to the late-time than the early Universe 
expansion history. Introducing the additional parameter 
weakens the constraints on most parameters consid- 



erably. This is also a consequence of the CMB data not 
being strongly sensitive to e(a). Therefore, other high- 
redshift probes such as 21 cm line spectra, gamma-ray 
bursts and Lyman-alpha forest data need to be added to 
improve the constraints for such models considerably. 

We next examine the dependence of our results on the 
form of the perturbations. It should be remarked that 
in Extended Quintessence (EQ) models, such as those 
represented by our 2-epsilons parametrization, the non- 
minimal coupling of the scalar field with the curvature 
leads to a specific behaviour of the perturbations, which 
is in general not equivalent to the one of a perfect fluid 
with unit rest-frame sound speed Cg = 1. This issue 
was discussed in great details in refs. [66] and [67]. In 
[66] it was found that, since the fluctuation in the en- 
ergy density of the scalar field, Sp-^^^ is directly related 
to the perturbation in the Ricci scalar, 5p-)^ tracks the 
density perturbation of the dominant component of the 
Universe. In MD, in particular, 5^ ~ 5m on all sub- 
horizon scales, (where the density contrast is defined by 
5% = 5 Pi/ 5 Pi) leading to the possibility of formation of 
dark energy clumps. In others words, in EQ Cs can drop 
much below unity, in contrast with what happens in nor- 
mal quintessence where = 1 such that dark energy 
perturbations are strongly suppressed on all sub-horizon 
scales. In principle, this can affect the growth of per- 
turbations in the matter sector and give observable ef- 
fects. In ref. [66], it was determined that values such as 

~ 0.001 during MD are typical for EQ models. More- 
over a small level of anisotropic stress is also present in 
EQ giving viscosity effects. In principle, the evolution of 
the perturbations in a given EQ model should be treated 
by solving the full system of perturbations equations for 
the specific choice of Lagrangian. Nevertheless, given the 
appropriate behaviour of the sound speed and the vis- 
cosity, the system can correctly be described in terms of 
fluid variables. To see how the constraints on our mod- 
els depend on the sound speed we carry out the analysis 




Figure 6: Posterior probabilities after marginalization over the other parameters: ACDM (dotted black), toCDM (dot-dashed 
blue), wZCDMl (dashed green) and ■u;ZCDM2 (solid red); all models with unit rest frame sound speed, Cs — 1. 



for several values for Cs, going as low as 10~^. On the 
other hand, as shown in [66], the anisotropic stress in 
EQ models is typically small, and so we neglect it. The 
results are shown in fig. 9, where we compare the pa- 
rameter estimation for = 1 with that for the extreme 
case of = 10^^. This quantifies the dependence of the 
results on our assumptions on the evolution of the DE 
perturbations. We see that for small Cg the value of ej\/ 
is less constrained, but in any case the bounds on eM 
and en remain of the same order of magnitude. The fact 
that the mean of the posterior of ej^ shifts to the right 
when allowing DE to cluster is in line with the tendency 
of the CMB data to favour an additional clustering rel- 
ativistic degree of freedom, see for instance [68]. The 
results for = 10"'^ are indistinguishable from those for 
= 10^^, so they are not shown. Thus, the model for 



the evolution of perturbations used in the previous anal- 
ysis is sufficiently general, at least at the present level 
of accuracy of the data. This is in agreement with the 
results of refs. [63, 69-71] who indeed show that, even 
for a small sound speed the effects are hardly detectable 
with current and even future data. 

IV. CONCLUSIONS 

In this paper we have studied a number of issues re- 
lated to vacuum fluctuations in a FRW background, and 
the role that they could play in explaining the observed 
acceleration of the Universe. 

First, we have studied vacuum fluctuations in a FRW 
background using a sharp cutoff in momentum space, and 
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Figure 7: luZCDMl: The marginalized likelihood contours at la and 2a are shown for all combinations of the parameters Us, 
QoE, Wo and e in toZCDMI. We compare ACDM (dotted black), wCDM (dot-dashed blue) and wZCDMl (dashed green) The 
shading is the mean likelihood surface for the wZCDMl model. All models have unit rest frame sound speed, Cs = 1. 



we have explicitly shown how one can recover the same 
results for the physical quantities that one obtains using 
a covariant regularization. The basic point is that, since 
a momentum-space cutoff breaks general covariance, one 
must allow for non-covariant counter-terms. A number 
of apparent pathologies discussed in the literature then 
disappear. In fact they are seen to be a consequence of a 
certain confusion between the bare and the renormalized 
(0|Tpi/|0). While the bare energy density and pressure 
do not satisfy energy-momentum conservation, we can 
renormalize the theory so that the renormalized quanti- 
ties do satisfy it. We have also clarified some confusion 
about the EOS parameter of the term proportional to 
H'^{t) in the vacuum energy density, for which general 
covariance requires w = Wtotit) rather than w ^ —1. 

Further, we have proposed an alternative to the stan- 
dard renormalization of vacuum energy, in which the 



A* divergence responsible for the cosmological constant 
problem is eliminated at the classical level, following a 
procedure that is already part of the standard tenets of 
classical general relativity and is used to define the ADM 
mass of a space-time. The crucial point here is to re- 
alize that, even if it is true that in GR the choice of 
the zero of the energy is not arbitrary, the correct choice 
is still non-trivial, and requires appropriate subtractions 
to be performed. We have seen that this point is fur- 
ther supported by the AdS/CFT correspondence, where 
the subtraction on the gravitational side corresponds to 
a renormalization of UV divergences on the CFT side. 
This subtraction scheme provides a possible solution to 
the "old" cosmological problem, namely why the vacuum 
energy density is not of order Mpj nor of order M**, where 
M is larger than the TeV scale, up to which quantum field 
theory is successfully applied. 
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Figure 8: toZCDM2: The marginalized likelihood contours at la and 2a are shown for all combinations of the parameters Us, 
f^DE, Too, Efl and e^f. Red solid lines are for u)ZCDM2 and blue dot-dashed lines are for wCDM. The shading is the mean 
likelihood surface for the U'ZCDM2 model. All models have unit rest frame sound speed, — 1. 



We then turned to the contribution to vacuum fluctua- 
tions due to an ultra-light scalar field with mass m < Hq. 
We have seen that, in general, such a field generates a 
contribution of the form f{x)R in the effective action 
(from which (0|T^^|0) is derived by functional differenti- 
ation). Such a term leads to cosmological models where 
the DE density tracks the dominant energy density com- 
ponent, thereby significantly alleviating the coincidence 
problem. We have performed a detailed comparison with 
cosmological observations of models where the DE den- 
sity has the form suggested by the above theoretical con- 
siderations. We have seen that, quite generally, we can 
describe these models with a two-parameter extension 
of wCDM cosmology and, using a full likelihood anal- 
ysis of CMB, BAO and SNIa data, we have presented 
the limits on the parameters and sm that encode the 
deviations from wCDM. They are consistent with their 
ACDM values, as it is also the case for wq in wCDM. 



Thereby there is at present no evidence for the presence 
of an ultralight scalar field that couples in this way to 
the curvature scalar. 
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Figure 9: Posterior probabilities after marginalization over the other parameters for ?i)ZCDM2 with: = 1 (soUd black) and 
Cs = 10~^ (dashed red). We note that for a small Cs, the value of €m is less constrained, a positive value is slightly preferred 
for ejj and that wq is less negative. 



Appendix A: Non-covariant counter-term for the 
H^{t)A'^ divergence 

Considering again for illustration a massless 
minimally-coupled scalar field, the bare energy den- 
sity and pressure, as far as the term cx H'^{t)A^ are 
concerned, can be written as 



Pbarc(Ac) 
Pbaro(Ac) 



167r2 
wtot(i) 



167r2 



(Al) 
(A2) 



where we used eq. (16). A counter-term proportional to 
the Einstein-Hilbert action. 



Sc^C{A,) I d^x^R, 



produces a counter-term 

2 ^5", 



'IT, 



jAU I 



(A3) 



(A4) 



where G^i, is the Einstein tensor. In FRW, Gqo = SH^ 



(where we used the relation 2H — — 3[1 + Wtot{'t)]H'^) 
This gives a counter-term for the energy density 



Pc = "6C(Ae)i?2(t), 
and a counter-term for the pressure 



(A5) 



(A6) 



and Gi 



The function C(Ac) can now be chosen so to cancel 
the divergence H^A^/{167t^) in eq. (Al) simply choos- 
ing 6C(Ac) = A2/(167r2) + ci, with ci a finite constant 
which corresponds to a finite renormalization of G, and 
can therefore be simply set to zero. This will also au- 
tomatically cancel the divergence WtotH^ {t)A'^/ {IQt^^) in 
eq. (A2). 

We then need a second counter-term which acts only on 
the pressure and cancels the term {2/3)H'^{t)A'^/{16TT'^) 
in eq. (A2). To explicitly write the required non- 
covariant counter-term, it is convenient to work directly 
in a FRW space-time with scale factor a{t). In order to 
perform the variation of the counter-term with respect to 
the FRW metric we introduce g^^i, ~ {—N'^{t),b'^{t)5ij) 
and perform the variation with respect to N{t) and b{t). 
Then, after performing the variation, we can set the lapse 
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function N{t) equal to one, and the function b{t) equal 
to the desired background value a{t). A counter-term 
'S'count in the action induces a counter-term for the en- 
ergy density given by 



Pcount 



2 ^'^count 

500- 



1 



while for the pressure 



SN 



N=l 



Pcount 



1 fss,. 



3a2 \^ Sb 



(A7) 



(A8) 



6(t)=a(t) 



The counter-term that cancels the divergence 
(2/3)iJ2(i)A2/(l67r2) in eq. (A2) will necessarily 
be non-covariant, and can in principle depend on the 
function a{t) both implicitly, due to the fact that the 



variation with respect to b{t) is evaluated in b{t) = a{t), 
and explicitly. A counter-term that does the job is 



967r2 
A2 
'l67r2 



dt —b\t) . 
a 



(A9) 



and iJ^iv and 



where = (-1,0^(5,^), = (-A''2 
R are the Ricci tensor and Ricci scalar constructed from 
5^1/ . By construction 5'count docs not depend on N(t), so 
it does not contribute to the energy density. Its contribu- 
tion to the pressure, from eq. (A8), is -AlH'^{t)/{24TT'^) 
and therefore cancels the term (2/3)if2(t)A2/(167r2) in 
eq. (A2). 
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